function [c,ceq] = fminconstr_irs(x,e2,ell,theta,sigma,pi,rho,eta,i,M,A,B)

c = []; % no nonlinear inequality
ceq = sys_irs(x,e2,ell,theta,sigma,pi,rho,eta,i,M,A,B); % the fsolve objective is fmincon constraints
